Finite-size Effects in a Two-Dimensional Electron Gas with Rashba Spin-Orbit 
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Within the Kubo formalism, we estimate the spin-Hall conductivity in a two-dimensional electron 
gas with Rashba spin-orbit interaction and study its variation as a function of disorder strength and 
system size. The numerical algorithm employed in the calculation is based on the direct numerical 
integration of the time-dependent Schrodinger equation in a spin-dependent variant of the particle 
source method. We find that the spin-precession length, L s controlled by the strength of the Rashba 
coupling, establishes the critical lengthscale that marks the significant reduction of the spin-Hall 
conductivity in bulk systems. In contrast, the electron mean free path, inversely proportional to the 
strength of disorder, appears to have only a minor effect. 
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I. INTRODUCTION 

The physical phenomenon behind the spin-Hall effect 
(SHE) 1,2 in two-dimensional (2D) systems is the flow of 
a pure spin current, spin-polarized in a transverse direc- 
tion, driven by a perpendicular electric field. Its exis- 
tence is conditioned by the presence of a spin-orbit in- 
teraction (SOI) , such as Rashba-Dresselhaus^ in n-type 
two dimensional systems or the spin-split band structure 
in p-type GaAs. 5 

If in clean samples the spin-Hall conductivity, <7 s h, 
was predicted to have a universal, constant value of 
e/87r, in the presence of disorder the resulting picture 
was less clear. It was pointed out that, in 2D infinite 
systems, in the presence of short range scatterers, the 
vertex corrections provided the exact compensation to 
cancel the effect^. Moreover, an argument was made 
that this cancellation occurs even for infinitesimal dis- 
order potentials 7 . These conclusions were challenged by 
analytic^ and numerical calculation a 9 ' 10 : 11 of the spin- 
Hall conductivity in one and two-dimensional finite-size 
mesoscopic samples, performed within the Landaucr- 
Biittiker formalism, where it was shown that the effect 
survives up to a critical disorder strength. 

Even though the robustness of the spin-Hall effect in 
the presence of disorder seems to have been definitively 
confirmed by the angle-resolved optical detection of spin 
polarization at opposite edges of a two dimensional hole 
layer—, a better understanding of the mechanism by 
which disorder and system size affect spin transport in 
systems with spin-orbit interaction warrants further in- 
vestigation. We focus therefore on a study of the in- 
terplay between the disorder strength, embodied in the 
electron mean free path I, and the spin precession length 
L s proportional to the spin-orbit interaction, in deter- 
mining the spin-transport regime in finite size samples. 
Such an analysis is especially relevant in two dimensions 
where, in the absence of any additional interactions, the 
two lengths are independent of each and along with the 



Fermi energy are the only relevant physical parameters 
of the system. 

The relationship between I and L s and the system size, 
L, determines the existence of four distinctive transport 
regimes. A semi-classical approximation is appropriate 
for L s <C L, when the spin coherence is lost over the 
length of the sample, while L s ^> L corresponds to a 
mesoscopic regime. When I 3> L the electron propaga- 
tion is ballistic, while for I <C L multiple scattering events 
are assumed and the diffusive regime is present. 

In the following analysis, we use the Kubo formula 
to estimate <t s h as a function of system size and disor- 
der in a two-dimensional electron system. The numerical 
formalism adopted here represents an extension to the 
spin-Hall problem of the particle-source method, devel- 
oped by Tanakaii This algorithm is based on the direct 
integration of the time dependent Schrodinger equation 
and allows the calculation of the matrix elements of the 
Green's functions, linear response functions, or any com- 
binations of Green's function and quantum operators in 
a very efficient way. 

The main result of this paper is that the delimitation 
between the mesoscopic and semiclassical regimes, as re- 
flected by the rapid decline of the spin-Hall conductivity, 
is established by L s . For system sizes smaller than L s 
the spin Hall conductivity increases monotonically with 
the system size, while being weakly affected by disor- 
der. When L 3> L s , <j s h decreases exponentially for any 
amount of disorder in the system. This result supports 
the conclusions of two previous reports by Sheng et alJ^ 
and Nomura et alJ^ where it was found that a S H remains 
finite up to an, unspecified, characteristic length scale 
and vanishes in the thermodynamic limit for any small 
amounts of disorder in the system. Here, we identify 
this length as being determined by the spin-precession 
length. Our results reflect no qualitative modification of 
the overall behavior when the system evolves from the 
ballistic to the diffusive regimes, crossover controlled by 
the mean free path characteristic lengthscale. For a fixed 
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Fermi energy and L s , the spin-Hall conductivity decrease 
monotonically with disorder for any system size, as the 
system evolves from diffusive to ballistic regime. 

II. THEORETICAL FRAMEWORK 
A. Model 

The single-particle Hamiltonian that describes the dy- 
namics of an electron of momentum p and effective mass 
m* , is written, in terms of the Pauli matrices a x . y and 
the Rashba coupling constant A, as: 

P 2 

H = — iA^-d^x). (1) 

The exact diagonalization procedure that can be per- 
formed on the Hamiltonian in the case of a clean system^, 
becomes impossible when disorder is included in the form 
of an additional random scattering term. It is therefore 
more convenient, for a numerical analysis, to adopt the 
tight-binding approximation for the many-body Hamilto- 
nian, by employing a local orbital basis associated with 
a virtual square N x N lattice, of constant a,Q. In this 
model, the many-body Hamiltonian, is: 

H = ^e iC \ a c ia -t J2 c L c j a (2) 

+ Vr [( c h c i+s*i - c !i c h-m) 

i,S x ,Sy 

In this expression, an electron with spin a at site i, cre- 
ated by c| Q , is subjected to a random on-site energy as 
in the Anderson model for disorder, generated by a box 
distribution £j e [—W/2,W/2\. The electron transport 
is described by a sequence of discrete hopping events. 
Lateral transport, without spin flip, to an adjacent site 
occurs with probability t = H 2 /2m* ao, taken to be the 
unit of energy in our calculation, as described by the 
second term in Eq. (J2J). Propagation along the diagonal 
sites, driven by the spin-orbit interaction, occurs with 
a simultaneous spin-flip, as in the last term of Eq. 
The latter is most important as it mixes the spin chan- 
nels and leads to a finite spin-Hall conductivity and spin 
accumulations at the edges of sample. The Rashba cou- 
pling constant is renormalized by the lattice constant to 
V R = HX/a Q . 

The Kubo formula for the spin-Hall conductivity is 
written as: 

asH = l Tr J^ ir^W) & [Gr{£) - Ga{£)] 

x v y G A (e) - f x G R (e) v y [G R (e) - G A {e)\) (3) 

The velocity operator is defined by the commutator: 
i hv y — [y, H], while for the spin current we adopt a tra- 



ditional expression^ given in terms of the anticommu- 
tator between the velocity operator and the Pauli ma- 
trix a z : j* = fi {a z , v x } /4. Gr/a(() represents the re- 
tarded/advanced Green's function. In Eq. © the inte- 
gration over the energy is restricted over the Fermi sur- 
face due to the presence of [— df{e)/de] factor. 

In the tight-binding framework, the effect of dis- 
order and spin-orbit interaction strength on the spin 
Hall conductance was investigated previously, using the 
Landauer-Biittiker formalis m 17 ' 18 . As will be discussed 
in the next section, in the present work we use a different 
approach for computing the Green's function needed for 
the calculation of the spin-Hall conductivity. 



B. Numerical algorithm 

Since the purpose of this investigation is an analy- 
sis of the spin-Hall conductivity dependence on system 
size and disorder, we will apply the Kubo formula to 
large size systems for different values of the disorder po- 
tential W. The numerical algorithm that underlies this 
calculation has been introduced in Ref. [13| and repre- 
sents an extension of the particle-source method com- 
bined with tight-binding formalism. This method was 
first applied to the calculation of the Green's function, 
density of states, conductivity^ and Hall conductivity^. 
The main advantage is that one can evaluate both the di- 
agonal and off-diagonal parts of the Green's function and 
their products with other quantum operators with a low 
computing effort. In principle the computing effort for 
computing the Green's function is 0(N 3 ), (Hamiltonian 
is expressed as a N x N matrix) while within the present 
algorithm only O(N) computational effort is required for 
the same calculation. Here, we briefly outline the main 
features of the algorithm. 

The central part of the method consists in solving 
the time dependent Schrodinger equation with a single- 
frequency source term: 

= H \l + li) *W ex?" , (4) 

at ' 

where rj is a finite small value and 9 is the step function. 
The solution of the equation, with the initial condition 
| j, t = 0) = becomes: 

\j,t) = -i f dt'e- iH ^\j)e- i ^ +i ^ t ' (5) 
Jo 

= 1 r-i(E+i n )t c - iHt )\j\ 

E + ir]-H l u/ ' 

For sufficiently large time, one can then write the solution 
to the Schrodinger equation in terms of the Green's func- 
tion acting on the "source" |j) with the relative accuracy 
6 = e-" T , as: 

\j,T) = G(E + ir,)\j)e-^ E +^ T , (6) 
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leading to a Green's function operating on the ket \j): 
G(E + ii]) \j) = lim \j, T) e l ^ E+l7l)T . (7) 

T — >oc 

The matrix element between states (i\ and \j) is then 
obtained as: 



(i\G(E + i V ) \j) 



lim 



It) 



(8) 



The matrix elements of a product including several 
Green's functions and other operators are obtained by 
choosing a new initial state, such as \f) = AG(E + irj) \j) 
in Eq. Q and repeating the same procedure. 

To calculate the matrix elements of the Green's func- 
tion at many different energy values, one solves Eq. Q 
simultaneously for a source term with multiple frequen- 
cies, \j) (J2i e~ l(El+i,l)t ) 6{t). Following the algorithm 
outlined above, one obtains as an approximate solution, 
the ket 



\J 



T)~^G(£ /+ «7) | j> e" 



■i(E,+iTi)T 



(9) 



The matrix element of the Green's function between the 
states (i\ and \j) for a giver energy is then easily obtained 
as 



TJo 



dt' (i\j,t')e i ^' +i ^ t ' , (10) 



where the terms involving transitions between different 
energies have been neglected with the relative accuracy 
S = 1/TAE, with AE the minimum increment in the 
energies E\. 

To obtain the time dependent ket \j;T), a direct nu- 
merical integration of the Schrodinger equation can be 
performed, as in the "leap-frog" algorithm 2 -. This is a 
second order, symmetrized differencing scheme, accurate 
up to (HAt) 2 . In this form, Eq. (|4]) becomes: 

\j;t + At) = -2iAtH\j;t) + \j;t- At) (11) 
- 2iAt\ ] )J2 e ~ l{El+u ' )t °( t ) > 



with a time step At determined by At = j3 / E max , where 
Emax is the absolute value of the extreme eigenvalue and 
(3 is a parameter whose value is less than 1 in order for 
the solution to be stable 2 ^. 

Estimating the trace in Eq. ((3J requires a suitable ba- 
sis set, such as the local orbital basis. It is more efficient 
however, to choose a randomized version of this basis, 

described by a ket \<f>) = Y^,n=i \ n ) ex P(~ * wnere 
\n) are the tight-binding orbitals and 4> n are random 
numbers in the [0, 27r] interval. For a given operator A, 
(<f>\A\)p) ~ (n| A \n) within the statistical errors of 

i/Vn. 



III. RESULTS AND DISCUSSION 

In this section we show the results of our computation 
based on the previously outlined algorithm. First, we 
study the Fermi energy dependence of the spin-Hall con- 
ductivity, presented in Fig. [TJ Random Fermi energies in 
the interval [-4t,4t] were considered and we average over 
2000 samples, for each system size. For clean systems 
and for states in the band, o s h is close to 0.8 (in unit of 
e/871"), except at half filling where due to electron-hole 
symmetry considerations it vanishes. In the presence 
of disorder, the calculated value of <j s h decreases, our 
results reproducing very well the known behavior previ- 
ously obtained in the Landauer-Biittiker formalism^iiS or 
by the analytical Kubo formula 15 . 
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FIG. 1: Spin Hall conductivity as function of Fermi Energy in 
units of e/8ir, for clean and disordered systems. System size is 
30x 30. Spin-orbit interaction strength is fixed to Vr = 0.2 1. 
Average was done over 2000 random frequencies. 



The dependence of the spin-Hall conductivity on the 
system size is shown in Fig. [5J One is interested 
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FIG. 2: Spin Hall conductivity as function of system size 
in units of e/8ir represented on a logarithmic scale. Fermi 
energy is fixed to Ef = —3.0 t and the spin precession length 
is L a ~ 16ao, corresponding to Vr = 0.2 1. Mean free path 
ranges from I ~ 45oo for W = 0.6 down to a few lattice 
constants when disorder increases. Averages are done over 
10 disorder samples and for each, 200 initial state vectors 
are considered. 
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in finding whether the variation of o s n is dramatically 
changed by disorder and what is the length scale at which 
this change occurs. For this, the two relevant param- 
eters are the electronic mean free path I and the spin- 
precession length L s . In a quasi-classical approximation 
I = 12fiVFao/{^N{EF) W) 2 , where vf is the Fermi ve- 
locity, while N(Ep) is the density of states at the Fermi 
energy measured from the bottom of the band. The spin- 
precession length is defined in terms of the Rasba cou- 
pling constant by L s = Tvtao/Vn. 

The electronic mean free path is the lengthscale that 
separates the ballistic from the diffusive regimes, with 
a ballistic behavior for system sizes smaller than / and 
diffusive otherwise. We found that the crossover be- 
tween these two regime is smooth, without any dramatic 
changes in the overall behavior of the spin-Hall conduc- 
tivity. The only observable effect is a decrease of the 
spin-Hall conductivity when disorder increases. For ex- 
ample, when W — 0.6 t, I ~ 45 do, while for W = 1.0 t, 
I ~ 16 ao, whereas, as can be seen in Fig. [5J the behav- 
ior of the spin-Hall conductivity remains unchanged. At 
the same time, for system sizes below L Sl o~ s h is always 
monotonically increasing, reaches a plateau between L s 
and 2L S , and then decreases for large system sizes, be- 
ing expected to vanish in the thermodynamic limit, as in 
Refs . 14 i 15 . The spin-precession length, therefore, is the 
characteristic lengthscale at which a crossover between 
the different regimes of the spin-Hall conductivity is ex- 
pected. In the semiclassical regime, a scaling analysis is 
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FIG. 3: Disorder dependence of the characteristic length scale 
£ for Fermi energy Ef = — 3.0 1 and spin-orbit interaction 
strength Vr = 0.2 1. 



appropriate. We find that, for a given Fermi energy, the 
size dependence of the spin-Hall conductivity can be very 
well fitted with a exponential function <j s h — exp(— 
where £ is a characteristic length that depends on the 
disorder strength, being divergent in a clean system. In 
Fig. [3] we present the dependence of this characteristic 
length as function of disorder which follows a power law 
£ ~ £ (t/W) - 57 behavior, with the best fit £ - 46 a . 
We remark again that the system size where <j s h starts 
to decrease is strongly conditioned by L s , rather than I as 
can be inferred from the similarities between the behav- 
iors drawn for different values of the disorder for the same 
L s . In all our calculations the decrease in the spin-Hall 
conductivity starts for system sizes L ~ 2L S regardless 
of the value of the electronic mean free path. 



IV. CONCLUSIONS 

In this work we study the effect of the spin-precession 
length scale and of the electronic mean free path on the 
the spin-Hall conductivity in different regimes, by adapt- 
ing the particle-source algorithm to spin transport in sys- 
tems with SOI in the framework of the the tight-binding 
approximation. The dependence of o~ s h on the Fermi en- 
ergy is also investigated. Our main finding is that the 
spin precession length is the critical length scale for the 
spin-Hall behavior. For a system size smaller that L s , 
the spin-Hall conductivity increases even in the presence 
of disorder, reaches a plateau between L s and 2L S and 
then in the semiclassical limit, when L L s decreases 
exponentially. In the thermodynamic limit, u s h is zero 
for any amount of disorder present in the system. We 
have also showed that the electronic mean free path does 
not play a fundamental role in the spin-Hall conductivity 
behavior. 
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